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P-f ABSTRACT 
O 

$_( ' Context. The weak gravitational lensing effect, small coherent distortions of galaxy images by means of a gravitational tidal 
field, can be used to study the relation between the matter and galaxy distribution. 

Aims. In this context, weak lensing has so far only been used for considering a second-order correlation function that relates 
• the matter density and galaxy number density as a function of separation. We implement two new, third-order correlation 
] functions that have recently been suggested in the literature, and apply them to the Red-Sequence Cluster Survey. As a step 
^ I towards exploiting these new correlators in the future, we demonstrate that it is possible, even with already existing data, to 

make significant measurements of third-order lensing correlations. 
, Methods. We develop an optimised computer code for the correlation functions. To test its reliability a set of tests involving 
■ mock shear catalogues are performed. The correlation functions are transformed to aperture statistics, which allow easy tests 
' for remaining systematics in the data. In order to further verify the robustness of our measurement, the signal is shown to 
J — . vanish when randomising the source ellipticities. Finally, the lensing signal is compared to crude predictions based on the 
' halo-model. 

J — ' Results. On angular scales between ~ 1' and ~ 11' a significant third-order correlation between two lens positions and one 
, source ellipticity is found. We discuss this correlation function as a novel tool to study the average matter environment of pairs 
of galaxies. Correlating two source ellipticities and one lens position yields a less significant but nevertheless detectable signal 
on a scale of ~ 4'. Both signals lie roughly within the range expected by theory which supports their cosmological origin. 
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1. Introduction 

The continuous deflection of light rays propagating 
through the Universe by inhomogeneities of the large-scale 
distribution of matter generates a coherent, weak distor- 
tion pattern over the sky: the shear field. As the source 



* Based on observations from the Canada-France-Hawaii 
Telescope, which is operated by the National Research Council 
of Canada, le Centre Nationale de la Recherche Scientifique and 
the University of Hawaii. 

Founded by merging of the Institut fur Astrophysik 
und Extraterrestrische Forschung, the Sternwarte, and the 
Radioastronomisches Institut der Universitat Bonn. 



of the shear field is gravity, the total matter distribution 
contributes to this effect. It is, therefore, a probe of the 
full matter content of the Universe. The shear field can 
be investigate d by correlating t he shapes of distant faint 
galaxies (e.g. Schneideil 2006t Van Waerbeke &: Mellier 



20031 : iBartelmann fc Schneideil l200lh. Since the first de- 
tection of this effec t (Bacon ct al. 2000: Kaiser et al.l l2000l : 
Van Waerbeke et al. 2000 ; Wittman et al. 2000. ) study- 
ing the weak gravitational lensing effect has become an 
important tool for cosmology. Together with the next- 
generation weak lensing surveys such as the RCS2 (on- 
going), CFHTLS (ongoing), Pan-STARRS 1 (commenc- 
ing soon) or KIDS (planned to commence during the 
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second half of 2008), the weak lensing effect will allow 
us to put tight con straints on cosmological parameters 
(|Peacocket al.ll2006l) . 

One important topic in contemporary cosmology is the 
relation between the dark matter and the galaxy popula- 
tion, the latter of which is thought to form under partic- 
ular conditions from the baryonic component within the 
dark matter density field. This relation can be studied by 
cross-correlating the shear signal and (angular) positions 
of a selected galaxy population. 

As the shear is quite a noisy observable, higher or- 
der galaxy-shear correlation functions are increasingly dif- 
ficult to measure. For this reason, studies in the past 
have focused on 2"'^-order statistics ( "galaxy-galaxy lens- 
ing", GGL hereafter) which involve one galaxy of the 
selected population (foreground) and one source galaxy 
(background) whose ellipticity carries the lensing signal. 
The GGL-signal can be used to learn more about the 
typical d ark matter environm e nt of si ngle galaxies (most 
rece ntly Kleinheinrich et al.' '20 06|; Man delbaum et al.^ 
2006b.a.G: .Seliak et al.. .2005. : .Hoekstra et al. 2005,. 2004c 
Sheldon et al.l 12004 ). or the so- called galaxy biasing 



2001|V. 



Simon et al.H2007t |Pen et al.ll2003l: [Hoekstra et aLll2002a. 



Schneider fc wattsi liooi) introduced "galaxy-galaxy- 



galaxy lensing" (GGGL) correlation functions and esti- 
mators thereof whi ch allow us to move to th e next, S'^'^- 
order level (see also I Watts fc Schneider 2005 ). The corre- 
lation functions now involve either two foreground galaxies 
and one background galaxy, or one foreground galaxy and 
two background galaxies. This idea was also discussed by 
I Johnstonl l 2006f l who studied how to derive the galaxy- 
galaxy-mass correlation function, which is one of the fore- 
going two, from weak gravitational lensing. These func- 
tions, although more difficult to measure than the two- 
point GGL signal, offer the opportunity to study the typ- 
ical environment of pairs of galaxies, e.g., within galaxy 
groups (or more technically, the occupation statistics of 
galaxies in dark matter halos, see .Cooray fc Sheth (200^ 
for a recent review), or possibly even the shape of dark 



matter haloes (ISm ith et al. 2006f l. More generally, they 



measure S'^'^-order moments between number densities of 
galaxies and the matter density of dark matter (cross- 
correlation bispectra). Hence, they "see" the lowest-order 
non- Gaussian features produced by cosmic structure for- 
mation. 

This paper applies for the first time the GGGL- 
correlation functions to existing data, the Red-Sequence 
Cluster Survey fRCS: iGladders fc Yeel[2005l ). and demon- 
strates that with the current generation weak lensing sur- 
veys it is already possible to extract these particular 3'''^- 
order statistics. 

The outline of the paper is as follows. We will give 
a brief description of the survey in Sect. [51 In Sect. [21 
we will define the correlation functions and their practical 
implementation as estimators for real data. In Sect. [4l 
our results will be presented, discussed and compared to 
halo-model based predictions to verify if the signal has 
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Fig. 1. Histogram of photometric redshifts of lenses, 
z G [0,0.4], and sources, z e [0.5,1.5], used in our anal- 
ysis. In total, we have got 2.4 x 10^ lenses {z « 0.30) and 
3.8 X 10^ sources (2 « 0.85). 

roughly the expected order of magnitude. Finally, in the 
same section, we demonstrate how the GGGL correlation 
function involving two lenses and one source can be used 
to map out the excess of matter - compared to the haloes 
of individual lenses - about pairs of lenses. 

Wherever a specific fiducial cosmology is needed 
^^m — 0.3, for the matter density parameter, and 
JIa — 0.7, for the dark energy density parameter, are as- 
sumed. Dark Energy is assumed to behave like a cosmo- 
logical constant. For the dark matter power spectrum nor- 
malisation we adopt erg = 0.9. 

2. Data: The Red-Sequence Cluster Survey 

The data used in this paper were t aken as part of 
the R ed-Sequence Cluster Survey (RCS; iGladders fc Yed 
I2OO5 ) , and comprise of approximately 34 square degrees of 
B, y, Rc and z' imaging data observed with the Canada- 
France-Hawaii Telescope (CFHT). The B and V bands 
were taken after completion of the original RCS, to allow 
for a better selection of clusters at low redshifts. These 
follow-up observations also enable the determination of 
photometric redshifts for a large sample of galaxies. This 
photometric redshift information is key for the work pre- 
sented here. A detailed discussion of these multicolour 
data, the reduction, and the photometric redshift deter- 
mination can be found in Hsie h et al. (2005). In the red- 
shift range out to z ^ 0.4 the photometric redshifts are 
well determined, with 70% of the galaxies within 0.06 of 
the spectroscopic redshift (as determined by comparing 
to a spectroscopic training set). For fainter galaxies the 
uncertainties become naturally larger. The photo-z uncer- 
tainty distribution in the RCSl photo-z catalogue is more 
or less a Gaussian for a given redshift range or a given 
apparent magnitude range. The relation between photo- 
z uncertainty, dz, and redshift is Sz ^ 0.06(1 + z). This 
relation over-estimates the uncertainty for z < 0.7 and 
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under- estimates it for z > 0.9 since the systematic error 
gets larger beyond that redshift. 

This photometric redshift catalogue was used by 



The galaxy-galaxy-shear correlator, 



Hoekstra et al.l (|2005D to study the virial masses of iso- 
lated galaxies as a function of luminosity and colour. To 
measure this galaxy-galaxy lensing signal, the photomet- 
ric redshift catalogue was matched against the catalogue 
of galaxies for which shapes were measured. This resulted 
in a sample of 8 x 10^ galaxies with 18 < Rg < 24 that are 
used in the analysis presented here. Hoekstra et al. ( 2005[ ) 
also present a number of lensing specific tests, demonstrat- 
ing the usefulness of the RCS photometric redshift cata- 
logue for galaxy-galaxy lensing studies. The frequency dis- 
tribution of photometric redshifts in our galaxy samples 
is shown in Fig. [TJ 

The galaxy shapes were determined from the Rc im- 
ages. The raw galaxy shapes are corrected for the ef- 
fects of the point-sprcad-function (PSF) as described in 
Hoekstra et al. (1998, 2002c). We refer the reader to these 
papers for a detailed discussion of the weak lensing anal- 
ysis. We note that the resulting object catalogues have 
been used for a range of weak lensing studies. Of these, 
the measurements of the le nsing signal ca u sed by large 
scale structure presented in iHoekstra et al. I (|2002bl l3l are 



particularly sensitive to residual systematics. The various 
tests described in these papers suggest that the systemat- 
ics are well under control. It is therefore safe to conclude 
that residual systematics in the galaxy shape measure- 
ments are not a significant source of error in the analysis 
presented here. 



3. Method 

Here we briefly summarise definitions of the three-point 
correlation functions, their estimators and the relation be- 
tween aperture statistics and correlation functions. A de- 
tailed derivation_aiid^ of those can be looked 
up in 



Schneider fc WattsI ^M)- 



3.1. GGGL-correlation functions 

For our analysis we consider two different classes of cor- 
relation functions. Both classes require triplets of galax- 
ies which are located at the positions 6i, 02 and 63 on 
the sky (see Fig. [2]). In a cosmological context, random 
fields - such as the projected number density of galaxies, 
N{9), or the shear field, 7(0) - are statistically homoge- 
neous and isotropic. For that reason, all conceivable cor- 
relations between the values of those fields depend merely 
on the separation, \6i — 6j\, and never on absolute posi- 
tions 6i. Therefore, our correlators are solely functions of 
the dimensions of the triangle formed by the galaxies. We 
parameterise the dimension of a triangle in terms of the 
lengths of two triangle edges, and 1^2, and one angle, 
03, that is subtended by the edges. Note that the sign of 
03, i.e. the handedness of the triangle, is important. 



0(^1,^2,03) = (Kg(ei)'^g(02)7 63; 



(1) 



is the expectation value of the shear at Os rotated in the 
direction of the line bisecting the angle 03 multiplied by 
the number density contrast of lens (foreground) galaxies 
at Oix- 



N(e) 

A rotation of shear is defined as 



7(0;^) 



(2) 



(3) 



where 7c is the shear relative to a Cartesian coordinate 
frame. It should be noted that Q and the following corre- 
lators are complex numbers. 

A second class of correlators are the galaxy- shear- shear 
correlators, 

G+(i?i,i?2,03) = (7(ei;¥'i)7*(e2;^2)'*g(e3)) , 
G_(7?i,z92,03) = {i{er,m)i{02;^2)K^{e3)) , (4) 

which correlate the shear at two points with the lens 
galaxy number density contrast at another point. Again, 
the shears are rotated, this time in the direction of the 
lines connecting the source (background) galaxies, at 61,2, 
and the lens galaxy at ^3. 

3.2. Practical estimators of correlators 

With practica l estimators for ([T]) and (|4]) in mind, 
Schneider fc W atts (200i) introduced modified correlation 
functions. They differ from Q and G± in that they are de- 
fined in terms of the number density of the lens galaxies, 
N{6), instead of the number density contrast. Kg: 



0(1^1,^2,03) 



{Nie^)Ni92h{9,;^)) 



(5) 



= 0(^9i,^2,03) + (7t) (^i)e-'^^ + (7t) (^2)e+'*^ , 
and 

G+(7?i,i92,03) = = (7(0i;¥'i)7*(^2;^2)A^(03)) (6) 



N 

G+(^?i, i?2, 03) + (7(^1; Vih*i02; V2)) 



G_(i?i, 792, 03) = = (7(^1; <Pi)7(02; ¥^2)^^(03)) (7) 
= G_(79i,792, 03) + (7(^1; ^1)7(^2; ^2)) . 

These correlators also contain, apart from the origi- 
nal purely 3'^'^-order contributions, contributions from 
2"'^-order correlations: (7t) {9) is the mean tangential 
shear about a single lens galaxy at separation 9 (GGL), 
(7(0i;(/3i)7*(02;</'2)) and (7(61; ^1)7(62; 952)) are shear- 
shear correlations which arc functions of the cosmic-shear 
correlators ^±{9) (e.g. Bartelmann & Schneider 20Ql|). To 
recover pure 3'''^-order statistics, the 2"'^-order terms can 
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Fig. 2. Geometry of the galaxy-shear-shear correlation, G+{ ^^,^2,(|)?.) {left pane l ), and the galaxy-galaxy-shear cor- 
relation, i92, 03) {right panel). The figure is copied from lSchneider fc Watts! ( 20051 ). 



either be subtracted, or even neglected, if we work in terms 
of the aperture statistics, as we will see in the next section. 

With respect to practical estimators, number densi- 
ties are more useful quantities because every single galaxy 
position is an unbiased estimator of N{9)/N. For that 
reason, every triangle of galaxies that can be found in a 
survey can be made an unbiased estimator of either Q 
(two lenses and one source) or G± (two sources and one 
lens). Since, generally, a weighted average of (unbiased) 
estimates is still an (unbiased) estimate^, we can combine 
the estimates of all triangles of the same dimension using 



arbitrary weights. 



for the sources. Note that for the 



following sums only triangles of the same "di, •§2 and 03 
have to be taken into account inside the sums. We adopt 
a binning such that "i^i, ^2 and 03 need to be within some 
binning interval to be included inside the sums, i.e. trian- 
gles of similar dimensions are used for the averaging: 



i,j.k—l 



J2 u'feeje^e-2"^^ e+2"^^ 



G 



1 

est _ hj,k=l 



(8) 



(9) 



where j,k £ {1...7Vs} are indices for sources and i G 
{1 . . . Ni} is the index of the lenses; A^i and Ns are the 
number of lenses and sources, respectively. By (pj and ipk 



^ The weighting scheme only influences the statistical uncer- 
tainty of the average, i.e. the variance of the combined esti- 
mate. Note that the whole statement requires that the weights 
are uncorrelated with the estimates that the average is taken 
of. 



we denote the phase angles of the two sources relative to 
the foreground galaxy i. 

The statistical weights are chosen to down-weight tri- 
an gles that contain sources who se complex ellipticities, 
Ci ( Bartelmann fc Schneidedliool . are only poorly deter- 
mined. Lenses, however, always have the same weight in 
our analysis. 

Similarly, we can define an estimator for Q. However, 
one has to take into account that ee"'^'^^"'"'^^-' of one single 
triangle - consisting of two lenses and one source with 
ellipticity e ~ is an estimator of 



{N{ei)N{92)-f{e3 



0) _ a(^i,^?2,03) 



{N{e,)N{e2)) 



1 + Lu{\e2 ~ Oil) 



(10) 



and not Q alone as has falsely been assumed in 
Schneider fc Wattd (l2005l) PI The function 



uj{\A0\) ^ {n^{e)n^{e + £^0)) (11) 

is the angular clustering of the lenses ( Peebles! 1980! ). 
Based on this notion, we can write down an estimator 
for^: 



E Wkeue-'^'P^+^'^\+u:{\e,-e,\ 

i,j,k=l 

(-1) E y^k 



(12) 



that includes explicitly the clustering of lenses. Here, Wk 
(fee {1 . . . iVs}) are the statistical weights of the sources. 
By ifi and ipj {i,j G {1 . . . TVi}) we denote the phase angles 
of the two lenses relative to the source fc. Again, only 
triangles of the same or similar dimensions (parameters in 
same bins) are to be included inside the sums. 



This becomes apparent if one sets 
7 {0^■, ieiiea) = 7 (6*3; ^g- ^) = const in the Eqs. (34) 
and (32), respectively, of ISchneider fc Watts! (|2005l ). 
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For obtaining an esti mate of coiO) in practice we em- 
ployed the estimator of Landv &: Szalavl ( 19931 ). which, 



transformations between correlation functions and power 
spectra. The aperture filter we employ for this paper is: 



compared to other estimators, minimises the variance to 
nearly Poissonian: 



, , DD DR 

uj{e) = 2 — 

RR RR 



1 



(13) 



It requires one to count the number of (lens) galaxy pairs 
with a separation between 9 and 6 + 59^ namely the num- 
ber of pairs in the data, denoted by DD^ the number of 
pairs in a random mock catalogue, RR, and the number 
of pairs that can be formed with one data galaxy and one 
mock data galaxy, DR. The random mock catalogue is 
computed by randomly placing the galaxies, taking into 
account the geometry of the data field, i.e. by avoiding 
masked-out regions, see Fig. [T^l We generate 25 random 
galaxy catalogues and average the pair counts obtained 
for DR and RR. 

When computing the Q and G± estimators, we suggest 
the use of complex numbers for the angular positions of 
galaxies: 'd = -di -\- 11^2 with di^2 being the a;/y-coordinates 
relative to some Cartesian reference frame (flat-sky ap- 
proximation). The phase factors turning up inside the 
sums (jH), (|9]) and (fT2|) are then simply (notation of Fig. 
El): 



-2i(^i 



where 'di 



13 . 

13 



23 . g-i(ipi + ip2) 
23 



1^13'^23 
1?13lll9 



231 



.(14) 



3.3. Conversion to aperture statistics 
In weak lensing, cosmological large-scale stru cture is often 



studi e d in terms of the apert ure st a tistics (ISimon et al 
20071: iKilbinger fc Schneided 120051 : IJarvis et al.l 12004 : 

ded fl 



Hoekstra et al.l l2002at ISchneideil (l998l: IVan Waerbekd 
1998 h that measure the convergence (projected matter 
distribution), k, and projected number density fields of 
galaxies. Kg, smoothed with a compensated filter u(x), i.e. 
/g dx xu{x) ~ 0: 



A4p(0) = ^^°°d2^u(^M^ K(|t9| 



(15) 
(16) 

where 9 is the smoothing radius. Map is called the aper- 
ture mass, while Af is the aperture number count of 
galaxies. With an appropriate filter these aperture mea- 
sures are only sensitive to a very narrow range of spa- 
tial Fourier modes so that they are extremely suitable 
for studying the scale-dependence of structure, or even 
the scale-dependence of re maining systematics in the data 
(jHetterscheidt et al.ll2007l ). Moreover, they provide a very 
localised measurement of power spectra (band power), in 
the case of (7V"M™) for n -I- m = 2, and bispectra, in 
the case of n + m = 3, without relying on complicated 



1 

2^ 



1 



X 

~2 



-2^2/2 



(17) 



as introduced bv ICrittenden et al. ( 20021 ). For an aper- 
ture radius of 9 the filter peaks at a spatial wavelength 
of ^ = which corresponds to a typical angular scale of 

27r TT /I 

As shear and convergence are both linear combinations 
of second derivatives of the deflection potential, the aper- 
ture mas s can be compute d from the shear in the following 
manner (jSchneider et al.lll998») : 



Map(0) + iAfi(f?) - ^ j^d^e'q(^^^^{e'-Ae'), (is) 
2 /"^ 

q{x) = — ds s u(s) — u(x) , (19) 
X Jo 

where we denote by ZO' the polar angle of the vector 6' . 
Note that in Eq. (fT8|) we place, for convenience, the origin 
of the coordinate system at the centre of the aperture. 



In expression (fT8| , Map is the E-mode, whereas M± is 
the B-mode of the aperture mass. Of central importance 
for our work is that we can extract E- and B-modes of the 
aperture statistics from the correlation functions. Since B- 
modes cannot be generated by weak gravitational lensing, 
a zero or small B-mode is an importan t check for a success 



ful P SF-correction of real data (e.g. iHetterscheidt et al 
20071). or the v iolation of parity-invariance in the data 



(ISchneidej[2003h . which is also a signature of systematics. 

Another argument in favour of using aperture statis- 
tics at this stage of our analysis is that 2"'^-order terms 
in Q and G+ do not contribute to t he 3'''*-order aperture 
statistics ( Schneider fc Wattsll2005l) . Therefore a signifi- 
cant signal in the aperture statistics means a true detec- 
tion of 3'^'^-order correlations. 

The 3'^^-order aperture statistics can be computed 
from Q via: 



(AA^Map) (01 



3) 



{N^Ma_){9i,( 



'2,03 



(20) 



(21) 



(l [g{di,d2,(^MNNM{^i,'92A: 



3|Cl, P2, 173 



where we have introduced for the sake of brevity an ab- 
breviation for the following integral: 



00 00 27r 

l[f]= J dM J d^2^2 J dh f 





(22) 



By 5R(a;) and 3(a;) we denote the real and imaginary 
part, respectively, of a complex number x. Eq. (|20[) is the 
E-mode of the aperture moment (A/'(^?l)A/'(^^2)Map(03)), 
whereas Eq. (|2ip is the corresponding parity mode that 
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is non-zero in the case of violation of parity-invariance; 
the latter has to be zero even if B-modes are present in 
the shear pattern that may be pro duced to some degree 
by intrinsic source alignment (e.g. Hevnians et al. 2004) 
or int rinsic ellipticity/shear correlations (jHirata fc Seliak 
|2004[) - that is, if we assume that the macroscopic 
world is parity- invariant. The integral kernel Aj^f^M for 
our aperture filter can be found in the Appendix of 
Schneider & WattsI (HqqI). 

The aperture statistics associated with the GGGL- 
correlator G± are the following: 



^[{MMAf) {01,62,93)^ 

(MlM) (91,62,03) - 
3?[(MM*A/-) (01,02,03) 



{MM*M) (01 



3)]/2, 



{MMU)(9i,02,03)]/2 



(M^M^pf^) (01,02,03) = 

9[(MMA/-) (01,02,03) + {MM*M) (01,02,03)]/'^ 
where we used the following definitions 
{MMM) (01,02,03) = 

I [g_(i9i, l?2, h)AMM^(^U^2, 4>3\01,02, 63)] , 

{MM*N) (01,02,03) = 

G+(^l,l92,h)AMM^M^l,d2, ^3\01,02, 0: 



(23) 



(24) 



(25) 



(26) 



(27) 



Eq. UMl) is the E-mode of {M,p(0i)M,p(e2)M(03)), Eq. 
((24|) is the B-mode which should vanish if the shear 
pattern is purely gravitational, and Eq. (|25p is again a 
parity-mode which is a unique indicator for systemat- 
ics. As before, the integral kernels Ama/tV and Amm'A/ 
for our aperture fil t er ma y be found in the Appendix of 
Schneider &: WattsI ()2005f ). 



3.4. Validating the code 

In the last section, we outlined the steps which have to 
be undertaken in order to estimate the 3"^-order aperture 
moments from a given catalogue of lenses and sources. 
The three steps are: 1) estimating the angular clustering 
of lenses yielding u!(9), 2) estimating Q and G± for some 
range of ■di^2 and for ip3 G [0, 2tt[, and finally 3) transform- 
ing the correlation function to (TV^Map) and (TVAffp) i^" 
eluding all E-, B- and parity- modes. 

There are several practical issues involved here. One 
issue is that, in theory, for the transformation we require 
g, G± for aU d e [0, oo], see Eq. (gl]). In reality, we will 
have both a lower limit (seeing, galaxy-galaxy overlap- 
ping), "iJiow, and an upper limit (finite fields), "iJuppor- On 
the other hand, the GGGL-correlators drop off quickly for 
large -d and the integral kernels Aj^j^m, Am mm, Amm'N 
have exponential cut-offs for di,d2 ^ ^1,2,3- Therefore, 
we can assume that there will be some range where we 
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12 5 10 

aperture radius [pixel] 

Fig. 3. Test run of our computer code with mock data 
based on some arbitrary convergence field. The mock 
data has been prepared such that {M^(0)AIap(0)) — 
{M(0)Mlp(0)) = {M^(0)); {NNN) = {j\f^e)) is the 
value that has to be found by the code (only equally sized 
apertures are correlated for test: 9i = 62 = 03 = 0). The 
binning range is between i? G [0.05, 200] pixel with 100 
bins; we use 10"* lenses and the same number of sources. 
For radii greater than ^ 2 pixel we get good agreement. 
The expected signal (solid line; computed from placing 
apertures) becomes inaccurate beyond 6' > 10 pixel be- 
cause the aperture size becomes comparable to the field 
size. The error bars denote the la sampling uncertainty 
due to finite galaxy numbers. The B- and parity modes 
(P) of the statistics are two orders of magnitude smaller 
than the E-modes and are oscillating about zero (plotted 
is modulus). 



can compute the aperture statistics with satisfactory ac- 
curacy. We perform the following test to verify that this 
is true: by using theoretical 3 D-bispectra of the galaxy - 
dark matter cross-correlations ( Watts &: Schneider 
we compute both the GGGL-correlation functions and 
the c orresponding aperture sta tistics (Eqs. 37, 38, 40, 51, 
52 of ISchneider fc WattsI l2005f ) . By binning the GGGL- 
correlators we perform the transformation including bin- 
ning and cut-offs in We find that one can obtain an 
accurate estimate of the aperture statistics within a few 
percent between roughly 0>4Oz?iow and < ??uppcr/10 
(using 100 log-bins for i?i_2 and 100 linear bins for (j)3). 
Therefore, with RCS-fields of typical size 139' we can ex- 
pect to get an accurate result between about 0'.5 < 6* < 
14'. 

Another issue is with step two above, in which the 
GGGL-correlators themselves need to be estimated. The 
estimators - Eqs. (jSj), ^ and (|12p - in terms of galaxy po- 
sitions and source ellipticities are simple but the enormous 
number of triangles that need to be considered is compu- 
tationally challenging (roughly lO^'^ per field for RCS). To 
optimise this process we employ a data structure based on 
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a binary tree, a so-c alled tree code (e.g. lJarvis et al.ll2004 : 
Zhang &: PenI 120051 ). The tree-code represents groups of 
galaxies within some distance to a particular triangle ver- 
tex as "single galaxies" with appropriate weight (and aver- 
age ellipticity). This strategy effectively reduces the num- 
ber of triangles. Moreover, we optimise the code such that 
only distinct triangles are found. Then, the other triangle 
obtained be exchanging the indices of either the two lenses 
(Q) or the two sources (G±) is automatically accounted 
for; this reduces the computation time by a factor of two. 

In order to test the performance and reliability of the 
code, we create a catalogue of mock data. In order to do 
this we use a simulated convergence field (K-field) on a 
grid, 512 x 512 pixeP, which has been obtained by ray- 
tracing through an N-body simulated universe. Actually, 
the only requirement that has to be met by the test field 
is that it behaves like a density contrast 6, i.e. {S) — 
and (5 > — 1, and that it has non- vanishing 3"^-order mo- 
ments, (6^) ^ 0. Based on this field wc simulate a shear 
and lens catalogue. The shear catalogue is generated by 
converting the K-field to a shear field and by randomly 
selecting positions within the field to be used as source 
positions. The positions and associat ed shear provide th e 
mock shear catalogue; for details see Simon et al.l (2004). 
In a second step, we use the K-field as density contrast, 
Kg, of the number density of lenses to make realisations 
of lens catalogues. This means one randomly draws posi- 
tions, 9, within the grid area and one accepts that posi- 
tion if X < ^+'^(^) where x is a random number between 
X G [0, 1] and K^ax the maximum value within the k field. 
Following this procedure one gets mock data for which 
K = Kg and therefore (7V"M™) = (A/'"+"'). In particular 
we must get, apart from the statistical noise due to finite 
galaxy numbers, (AA2(0i)Map(02)) = {J^i0i)Mipi02)) 
when running our codes with the mock data. 

Parallelly, we smooth the test shear field within aper- 
tures according to the definitions (|18p with our aperture 
filter and estimate the test data aperture statistics directly 
by cross-correlating the smoothed fields. This also has to 
be comparable (apart from shot noise) to our code output. 
The result of this test can be found in Fig. [31 

As a further test we take the same mock data but 
rotate the ellipticities of the sources by 45 degrees, i.e. 
we multiply the complex ellipticities by the phase factor 
_g-2i0 -^^ith <j) — 45°. This generates a purely B-mode sig- 
nal that should only be picked up by the B-mode channels 
of the aperture statistics, yielding a plot similar to Fig. [3] 
The parity mode in (A/'M^p) has to be unaffected. This is 
indeed the case (figure not shown). 

The test results make us confident that the computer 
code is working and that we achieve a good accuracy even 
though we are forced to bin the correlation functions and 
to use a tree-code that necessarily makes some additional 
approximations . 
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Fig. 5. Combined measurement of angular clustering 
of our sample of lenses (no correction for the integral 
constraint). Error bars were obtained by looking at the 
field-to-field variance. The solid line is a power-law fit, 
uj{0) = A^e-I^, to the regime 9 G [0'.1,3']. 

4. Results and discussion 

We applied the previously outlined method to the RCS 
shear and lens catalogues. 

Lenses were selected between photometric redshifts 
< z < 0.4, whereas sourc es were from t he r ange 
0.5 < z < 1.4. Compared to iHoekstra et al.l (j2005l ). in 
which photometric redshifts smaller than 0.2 were ex- 
cluded, we were less strict about the lowest redshift of the 
lenses. This is likely to have introduced some misidentified 
lenses into our sample (less than 10%, see Fig.[T]) as RCS is 
lacking a U-band filter. Moreover, including sources with 
photometric redshifts larger than z ~ 1.0 is also rather op- 
timistic because photometric redshifts within that range 
can become quite unreliable as well. Therefore the tail of 
the redshift distribution in Fig. [T] may be slightly inac- 
curate. Still, sources with photo-z's greater than one are 
likely to be high-redshift galaxies. 

However, for our purpose, namely demonstrating a ro- 
bust detection of GGGL, the biases in the redshift distri- 
bution of lenses and sources are acceptable. These biases 
in the estimated redshift distribution only become an issue 
if one wants to thoroughly model the GGGL-signal. 

4.1. Aperture statistics 

As a first result we would like to draw the reader's at- 
tention to the angular clustering of lenses which is plot- 
ted in Fig. [5l This measurement was required for the es- 
timator Q in Eq. (fT2|) . As widely accepted, the angular 
correlation function w(0) is, for the separations we are 
considering here, well approximated by a simple power- 
law, depending on gal axy type, colour and luminosity (e.g. 
Madgwick et al.|[20 03V As can be seen in Fig.[5l the power- 
law behaviour is also found for our lens galaxy sample. The 
angular clustering plotted is still affected by the so-called 
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Fig. 4. Left: Aperture statistics (M'^M^.p) {0, 6, 9) for different aperture radii 6 as measured in RCS. The upper panel 
is the E-mode, the lower panel is the parity mode which is consistent with zero. Error bars denote the field-to-field 
variance between the ten RCS fields. Statistical errors are strongly correlated. The lines are tentative halo model-based 
predictions with arbitrary HODs for a ACDM cosmological model (see text). Right: Aperture statistics (A/'M^p) (0, 9, 9) 
for different aperture radii 9 as measured in RCS. The upper panel contains the E-mode measurement, while B-mode 
(stars) and parity mode (squares) are plotted inside the lower panel. Error bars that extend to the bottom of the 
upper panel denote data points that are consistent with zero. 



integral constraint ( Groth h Peebles! 1X977 ). which shifts 
the estimate of w downwards by a constant value depend- 
ing on the geometry and size of the fields. For small ^ ^ 3' 
this bias is negligible so that we used only the regime 
9 S [0^1, 3'] to find the maximum likelihood parameters of 
the power-law. 

For Q this power-law fit was used. Possible deviations 
of the true clustering from a power-law for 9 > 2' were 
negligible because for the estimator one actually needs 1 -f 
u) instead of lo. Since lu is roughly smaller than ^ 0.05 and 
decreasing for 9 > 2' , we gather that a certain remaining 
inaccuracy in to has no big impact on 1+uj. The power-law 
index is, with /3 — 0.58, fairly shallow, w hich is typical for 



a rela tively blue sample of galaxies (e.g. iMadgwick et al. 
2003h . 



In a second step, the correlation functions Q and G± 
were computed separately for each of the ten RCS fields. 
The total combined signal was computed by taking the 
average of all fields, each bin weighted by the number of 
triangles it contained. For the binning we used a range 
of 0".8 <-& <54' with 100 bins, thus overall 10^ triangle 
configurations. By repeatedly drawing ten fields at ran- 
dom from the ten available, i.e. with replacement, and 
combining their signal we obtained a bootstrap sample of 
measurements. The variance among the bootstrapped sig- 
nals was used to estimate the sum of cosmic variance and 
shot noise, thus the remaining statistical uncertainty of 
the correlation functions. 

Finally, the correlation functions were transformed to 
the aperture statistics considering only equally sized aper- 
tures, i.e. 6*1 =02 = ^'3, see Fig. |4l For the scope of this 
work, equally sized apertures are absolutely sufficient. In 



future work, however, one would like to harvest the full 
information that is contained in these statistics by explor- 
ing different 9i which then would cover the full (projected) 
bispectrum. 

For a start, we would like to focus on Q. The left panel 
in Fig. d] reveals a clean detection of ^A/'^Afap) for aper- 
ture radii between 0'.5 < 9 < 5' (with the adopted filter 
this corresponds to typical angular scales between 1' and 
11') demonstrating the presence of pure 3'^'^-order corre- 
lations between shear and lens distribution in RCS. The 
parity mode of this statistic is consistent with the zero as 
expected. Fig. |4] is one of the central results of this paper. 

We would like to further support that this is a real, 
i.e. cosmological, signal by comparing th e measurement to 
crude halo model-based predictions (see ICoorav &: Sheth 
2OO2I for a review). The halo model was used to pre- 
dict a spatial cross-correl ation bispectrum, Bggs (Eq. 12 in 



Schneider fc mttiliooi), for a particular fiducial cosmo- 



logical mode l and halo occupation dist ribution (HOD) of 
galaxies (see lBerlind fc Wei nberg 2G0g) By applying Eqs. 
(21), (52) in ISchneider fc Watts (,2005h . Bggs was trans- 
formed, taking into account the correct redshift distribu- 
tion of lenses and sources (Fig. [1]), to yield the aperture 
statistic s. A standard concor dance ACDM model was em- 
ployed (jBardeen et al. Il98d) with parameters J^a — 0.7, 
for the dark energy density, Om = 0.30, for the (cold) dark 
matter density, erg = 0.9 for the power spectrum normal- 
isation, and r — 0.21 for the shape parameter. This is in 
agree ment with constrain ts based on the first WMAP re- 



lease (jSpergel et al.ll2003l) 



The latest constraint s favour a somewhat smaller value 
for fTg (jBeniamin et al 



20071 : iHetterscheidt et al.l l2007^ 
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Fig. 6. Residual signal (squares) of GGGL in RCS when the ellipticities of the sources are randomised {left: (N'^M^p^, 
right: I^MM'^^). For comparison, the original signal before randomisation is also plotted (crosses). The line is a crude 
halo-model prediction of a blue galaxy population as in Fig. |4l The error bars of the randomised signal quantify the 
background noise of a null-signal. This indicates that we have a significant detection of (A/'^Map) in the left panel but 
only a weak detection of (A/'M^p), most significant at about 2', in the right panel. 



which would shift the expected amplitude of GGGL to- 
wards smaller values . If we apply the scaling relation of 
Jain fc SeliakI (1997), given for the convergence bispec- 
trum, as a rough estimate of this shift, B^, oc CTg'^, we ob- 
tain a correction factor of about two for erg = 0.8 {N and 
Map should have the same erg-dependence for unbiased 
galaxies) . 

The halo-model predictions depend strongly on the 
adopted H OD. The basic set up for this model was that 
outlined in' Takada fc Jain ( 2003*) . which splits the occupa- 



tion function, N{M), into contributions from "red", TVfl, 
and "blue" ,Nb, galaxies: 



{Nb) [M] = 



TUB 



IB 



-I- Aexp (--Ao(logio("i) -tobJ^) 



(Nr) (M) = — 



exp 



1,2 



(28) 



As parameters we used niB ~ 2.34 x 10 M©, A = 
0.65, Ao = 6.6, niB, = 11.73, = 1.8 x IO^^Mq and 
TOi^Q = 4.9 X 10^^ Mq. Blue galaxies have a peak halo 
occupancy of around lO'^^ Mq and a shallow power law 
(7s = 0.93) at high halo masses. In this simple prescrip- 
tion, red galaxies are relatively more numerous in higher 
mass halos (7^ = 1.1) and are excluded from low mass ha- 
los by an exponential cutoff around 5 x 10^^ M©. Factorial 
moments of the occupation distribution - the cross bis- 
prectra Bggi, and Bddg require the mean and variance 
- were as prescribed in the model of iScoccimarro et all 
( 20011 ). In this way, the moments are Poissonian for 
higher mass halos, becoming sub-Poissonian for masses 
below 10^3 Mr:,, i.e. (N'^) [M) = oP-{{N) {M)f , where 



^0 

a = O.Slogio (m/lO^^M©). 

We stress at this point that we made no attempt 
to "fit" parameters to the data, we merely intended to 



bracket a range of possible results. To choose a range of 
plausible scenarios, we constructed the theoretical aper- 
ture statistics for "red" galaxies, "blue" galaxies and for 
"all" galaxies (in which the occupation functions for red 
and blue galaxies are added together directly). We also 
showed predictions for the unbiased case, in which the oc- 
cupation function N{M) oc M with Poisson moments for 
<^7V^(M)). Galaxies were assumed to follow the CDM halo 
density profile (NFW) with no assumption of a central 
galaxy. Other parameters that define the halo model set 
u p (e.g. concentration o f the NFW profile) were as used 



Takada fc Jah] (|2003l) . 



Our measurement of (AA^Map) lies somewhat above 
the lower bound of the expected physical range of values, 
giving support as to the cosmological origin of the signal. 
Moreover, taken at face value, our result appears to fit 
the picture that the lens population consists of rather blue 
galaxies as has been concluded from the shallow slope of 
the angular correlation function uj. 

We randomised the ellipticities of the sources and re- 
peated the analysis. Since the coherent pattern, and its 
correlation to the lens distribution, is responsible for the 
signal, destroying the coherence by randomising the ellip- 
ticity phase should diminish the signal. That this is the 
case can be seen in Fig. [HI (left panel). 

Analogous to (TV^Map) we computed and predicted 
(A/'Mfp), the result for which is shown in Fig. 2] (right 
panel). Here a signal significantly different from zero was 
only found for aperture radii 1' < < 3' and at about 
9 ~ Below 6 ^ I'.b the parity mode is not fully consis- 
tent with zero. Hence, we may have a non-negligible con- 
tamination by systematics in the PSF correction and/or 
intrinsic alignments of the sources that may hamper a 
clean detection. For radii where we find a non-zero signal. 
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Fig. 8. Signal in Q originating from pure 2"'^-order statis- 
tics (GGL) that was subtracted from Q to obtain Fig. [9] 
(left). The gray-scale intensity is the modulus of the 
sticks indicate the average shear at the source position 
in the lens-lens-source triangle configuration. The units of 
the axis are in /i~^Mpc which corresponds to the mean 
physical scale at the lens distance of about z = 0.30. The 
two lenses are located at the positions of the crosses, left 
and right from the centre. 



the signal is on average smaller than the lowest theoretical 
value from our crude models. However, as discussed above, 
a lower erg easily brings the model down towards smaller 
values. The signal disappeared if the cUipticitics of the 
sources were randomised (Fig. [6l right panel). Therefore, 
we found a tentative detection of (A/'Affp) in our data. 



4.2. Mapping the excess matter distribution about two 
lenses 

The aperture statistics clearly have advantages: the B- and 
parity-modes allow a check for remaining systematics in 
the data, and 2'^'^-order statistics do not make any contri- 
butions so that we can be sure to pick up a signal solely 
from connected 3'^'^-order terms. This is what we did in 
the forgoing subsection. The result suggests that we have 
a significant detection of G- 

The disadvantage of using aperture statistics is, how- 
ever, that they are hard to visualise in terms of a typi- 
cal (projected) matter distribution (lensing convergence) 
about two lenses, say. Therefore, we introduce here an al- 
ternative way of depicting G which is simil ar to the work 
that has been proposed by Johnston! ( 2006h . 

A similar way of visualising G± probably could be 
thought up as well. However, since we found only a weak 
detection of GGGL with two sources and one lens we post- 
pone this task to a future paper and focus here on Q alone. 

The following summarises what essentially is done if we 
estimate Q from the data for fixed lens- lens separations. 
We pick out only lens-lens-source triangles from our data 
set in which the lenses have a fixed separation or a sepa- 
ration from a small range. Each triangle is placed inside 
the plot such that the line connecting the lenses is paral- 
lel to the z-axis and that the centre of this line coincides 
with the centre of the plot, as seen for the triangles in Fig. 
[71 The ellipticities of the sources of all triangles are then 
multiplied by 1 -I- wd 6*2 — 6*1 1) (rescaled according to Eq. 
[TUl) and (weighted) averaged at the source positions. For 
this paper, we used 128 x 128 grid cells for binning the el- 
lipticities. Following this procedure we effectively stacked 
all shear patterns about a lens-lens configuration - rotated 
appropriately - to obtain an average shear field about two 
lenses. This is, in essence, the meaning of Q. The full ^ is a 
bundle of such plots with continuously changing lens-lens 
separations. 

Note that the ellipticity at the source position, stored 
in Q, is rotated by <j>3/2 (Fig. [21 right panel). For the 
following plots, on the other hand, we used the shear in 
Cartesian coor dinates relative t o axis defined by the lens 
positions, as in I Johnston ( 20061 ). Therefore, when gener- 
ating the plot we were rotating our measurements for Q 
appropriately. 

The resulting plot has symmetries. Firstly, we do not 
distinguish between "lens 1" and "lens 2" . Both lenses are 
drawn from the same galaxy sample. This means for every 
triangle, we will find the same triangle but with the posi- 
tions of "lens 1" and "lens 2" exchanged. Therefore, the 
two lenses and the source of the triangle named "original" 
in Fig.[7lwill make the same contribution but complex con- 
jugated at the source position of the triangle named "in- 
dex swapping" . Thus, quadrants I and III will be identical 
apart from a complex conjugate and mirroring the posi- 
tions about the x- and y-axis. The same holds for quad- 
rants II and IV. This would no longer be true, of course, if 
we chose the two lenses from different catalogues in order 
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Fig. 9. Plots of Q after subtraction of the 2"'^-order signal from Q. The units used are /i"^Mpc, which corresponds 
to the mean comoving physical distance at the lenses' distance of, on average, z = 0.30. Le/i; Lenses were selected to 
have a mutual angular separation between 40" and 80" corresponding a projected physical scale of about 250 /i~^kpc. 
Right: Lenses were chosen to have a separation between 4' and 8', or equivalently a projected comoving separation 
between 1 — 2 h^^Mpc. 



to, for instance, study the matter distribution around a 
blue and a red galaxy. 

A second symmetry can be observed if the Universe 
(or the PSF-corrected shear catalogue) is parity invariant. 
Mirroring the triangle "original" with respect to the line 
connecting the two lenses (a;-axis) results in another tri- 
angle coined "parity" . For parity invariance being true the 
ellipticity at the source position of "parity" is on average 
identical to the ellipticity at the source position of tri- 
angle "original". In this case, quadrant IV is statistically 
consistent with quadrant I and quadrant II with quad- 
rant III (after mirroring about the z-axis). Taking parity 
symmetry for granted could be used to increase the signal- 
to-noise in the plots by taking the mean of quadrants IV 
and I (or II and III). 

Since the way of binning in the plot is completely dif- 
ferent from the way used to get the aperture statistics 
out of RCS, we made two reruns of the estimation of G 
with our data. For the first run wc only considered lens- 
lens separations between 40" and 80", the second run se- 
lected triangles in which the lenses had a separation be- 
tween 4' and 8'. For a mean lens redshift of z ^ 0.3 this 
corresponds to a projected physical comoving separation 
of roughly 250/i~^kpc and 1.5 ft.~^Mpc, respectively. As 
usual, the results from the ten individual fields were av- 
eraged by weighting with the number of triangles inside 
each bin and the statistical weights of the sources. 

Since we effectively stacked the shear fields about all 
pairs of lenses, aligned along the lens-lens axis, we ob- 
tained the average shear about two lenses. The shear pat- 
tern still contained a contribution stemming from GGL 
alone. This contribution could, however, easily be sub- 



e.g. 

2nd. 



tracted according to Eq. ^ after estimating the mean 
tan gential shear (7t) W , about single lens galaxies (see 
Simon et al.l l2007l) . A typical shear pattern due to 
order GGL can be seen in Fig. [51 This is the shear 
pattern that is to be expected if the average shear about 
two lenses is just the sum of two mean shear patterns 
about individual lenses. They contain all contributions 
that are statistically independent of the presence of the 
other lens. Therefore, contributions (contaminations) to 
the shear from lens pairs that are just accidentally close 
to each other by projection effects, but actually too sepa- 
rated in space to be physically connected, are removed. 

Now, Fig. [5] shows the shear patterns after removing 
this signal. Clearly, there is a residual coherent pattern 
which is most pronounced for the smaller lens-lens sep- 
arations. This proves that one finds an additional shear 
signal around two galaxies if they get close to each other. 
Hence, the average gravitational potential about two close 
lenses is not just the sum of two average potentials about 
individual lenses. 

Unfortunately, all physically close galaxies with a fixed 
projected angular separation contribute to the excess 
shear- independent of whether they are in galaxy groups 
or clusters. Exploiting lens redshifts and rejecting lenses 
from regions of high number densities on the sky might 
help to focus on galaxy groups, for example. This, how- 
ever, is beyond the scope of this paper. 

One can relate the residual shear pattern in Fig.[9]to an 
excess in projected convergence (matter density) using the 
well known relation between convergence and cosmic shear 



12 



Simon et al.: GGGL in RCS 




Fig. 10. Convergence fields obtained by transforming the shear fields in Fig. [9l They are related to the (average) 
excess in matter density around two galaxies of fixed angular separation after subtraction of the matter density profile 
that is observed about individual galaxies. Left: Residual convergence for two lenses with projected comoving distance 
of roughly 250/i~^kpc. The box-size is 1.7ft.~^Mpc x 1.7 /i~^Mpc. Right: Residual convergence at about 1.5/i^^Mpc 
projected lens-lens distance. The box-size is 8.5ft.~^Mpc x 8.5 h~^Mpc. Note that the convergence in this figure is 
lower by roughly an order of magnitude compared to the left figure. 



in weak gravitational lensin g (jBartelmann fc Schneider 
200lUKaiser fc Sauiresl[l993h : 



(29) 



where 7^ and Ki are the Fourier coefficients of the 
shear and convergence fields, respectively, on a grid and 
^ = (£1,^2) is a particular angular mode of the grid in 
Cartesian coordinates. We obtained the 'fg's by employ- 
ing Fast-Fourier- Transforms (and zero-padding to reduce 
undesired edge effects) after binning the residual shear 
patterns onto a 512 x 512 grid. We assumed that the con- 
vergence is zero averaged over the box area which makes 
Ke = for i = 0. Fig. \TU\ shows the thereby computed 
maps. The plots were smoothed with a kernel of a size 
of a few pixels. 

As a cross-check we also transformed the shear pat- 
tern produced by the 2"'^-order terms in Q (Fig. [8]) and 
found, as expected, that the corresponding convergence 
fields were just two identical radially symmetric "matter 
haloes" placed at the lens positions in the plot. 

In the same way as in the previously discussed shear 
plots, parity invariance can also be checked in the conver- 
gence plots: quadrants I and IV (or II and III), mirrored 
about the x-axis, have to be statistically consistent. If we 
would like to enforce parity invariance, we could take the 
average of the two quadrants. Secondly, if one obtains the 



convergence field from the shear field via a Fourier trans- 
formation as described before, the convergence field will 
be a field of complex numbers. In the absence of any B- 
modes, however, the imaginary part will be zero or pure 
noise. Thus, the imaginary part of the convergence can be 
used to either check for residual B-modes or to estimate 
the noise level of the E-mode (real part). 

This was done for Fig. [TlJ We found that the residual 
convergence for the small lens-lens separation is highly 
significant within the central region of Fig. [TOl left panel, 
whereas the convergence in the right panel of Fig. \W\ is 
noise dominated. This means we did not find any excess 
convergence beyond the noise level for the lens-lens pairs 
of large separation. 

To sum up, one can see that closer lens pairs are 
embedded inside a common halo of excess matter, while 
the lenses with larger separation appear relatively discon- 
nected; the convergence for the lenses of larger separation 
is lower by at least one order of magnitude and slips be- 
low noise level in our measurement. This result definitely 
deserves further investigation which we will do in a forth- 
coming paper. 

5. Conclusions 

We found a significant signal of GGGL in RCS - at least 
for the case for which we considered two lenses and one 
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Fig. 11. Plots similar to the plots in Fig. [TOl except the the shear has been rotated by 45° (B-mode) before transforming 
to the convergence fields. The thereby obtained convergence quantifies the statistical noise in the plots of Fig. [TOl 



source. The signal is of an order of magnitude which is 
expected from a crude halo model-based prescription. This 
suggests a cosmological origin of the observed correlation. 

In particular, our finding demonstrates that wide-field 
surveys of at least the size of RCS allow us to exploit 
GGGL. As can be seen in Fig. |4] (left), the remaining sta- 
tistical uncertainties of the measurement are much smaller 
than the shift of the signal expected for different HODs of 
the adopted halo model. This means with GGGL we now 
have a new tool to strongly constrain galaxy HODs, and 
possibly even spatial distributions of galaxies inside haloes 
in general, which is a parameter in the framework of the 
halo model. As the wide-field shear surveys of the next 
generation will be substantially larger than RCS those 
constraints will become tighter. Further subdivisions of 
lens samples into different galaxy types and redshifts will 
therefore still give a reasonable signal-to-noise ratio. 

Leaving the interpretation in the context of the halo 
model aside, the measurement of GGGL can be trans- 
lated into a map of excess convergence around two galax- 
ies of a certain mutual (projected) distance. For RCS, we 
demonstrated that there is a significant excess in conver- 
gence about two lenses if galaxies are as close as roughly 
250/i^^kpc. Although the details need still to be worked 
out, this promises to be a novel way of studying the matter 
environment of groups of galaxies. 
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Fig. 12. Determined geometry of the individual RCS fields (field masks). Each mask has a size of approximately 
139' X 139' or, cqiiivalently, 20.000 x 20.000 pixeP. Black colours denote masked-out regions, while white indicates 
regions where galaxies have been found. 



